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SUMMARY 



Current methods that take into account interference 
and compressibility effects have been used to estimate 
pressures over and near turrets located on a wing in 
regions of high interference velocities. 

Graphical comparisons of measured pressures, avail- 
able from wind-tunnel tests, and estimated pressures are 
presented for angles of attack of 0°, 1°, and 4° at Mach 
numbers up to 0.650. The agreement between measured and 
estimated pressures indicated that present methods are 
adequate for predicting approximate critical speeds and 
for estimating loads at subcritical speeds with suffi- 
cient accuracy for design purposes. The methods used in 
this paper were not applicable at speeds much higher than 
the critical. 



INTRODUCTION 



The investigation reported in the present paper 
was made in order- to check the applicability of current 
methods, suggested in reference 1, to pressure estima- 
tions involving large interference effects. Prior to 
this investigation, successful predictions of pressure 
distributions for gun turrets and similar shapes were 
for the most part limited to protuberances located on 
fuselages, where the interference was usually negligible. 
The applicability of current methods to estimations 
involving high interference velocities was uncertain 
because of insufficient experimental data with which to 
compare estimates. Recent tests of a turret model in 
the NACA 3-foot high-speed tunnel (unpublished), 
however, provided the experimental data necessary for 
investigating the accuracy of predictions for turret-like 
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protuberances located in regions of high interference 
velocities. The purpose of the present paper is to 
apply current methods, including interference and 
compressibility effects, to the estimation of pressure 
distributions for this turret, for which measured data 
are available, and to compare the estimates with meas- 
urements. 

Present theory, which does not hold at speeds much 
higher than the critical, limited the Mach number range 
for the estimations. Practical comparisons were made at 
a subcritical speed (Mach number M ~ 0.275) and at the 
approximate critical speed (M =0.500); a few comparisons 
were made at a supercritical speed (M = 0.650) to show 
the inadequacy of the nethods used at subcritical speeds 
when applied in the supercritical range. Angles of 
attack for which the estimations were made (0°, 1°, and 
4°) were determined largely by the range of experimental 
values available for comparison. 

This and similar applications of current esti- 
mating methods are valuable because they indicate that 
predictions may be substituted for the more expensive 
wind-tunnel tests of turret-like protuberances when 
approximate critical speeds and loads at subcritical 
speeds are desired. 



SYMBOLS 



P 



pressure coefficient 



V 



velocity 



v/v 0 



velocity ratio; total unless indicated otherwise 
by subscript 



velocity increment 



AV/V 0 



velocity- increment ratio, induced velocity, or 
interference velocity 



M 



Mach number 



a 



wing angle of attack,, degrees 
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Y ratio of specific heat at constant pressure to 

■ specific heat at constant volume for air 

( c p/ c v) 



Sub scripts? 


aos 


absolute 


cr 


critical; at speed of sound 


i 


inc. omp r e s s ib le 


max 


maximum value 


min 


minimum value 


o 


free stream 


t 


turret 


w 


wing 


2D 


two-dimensional flow 


3D 


three-dimensional flow 



DESCRIPTION OF TURRET SHAPES 



The turret installation, for which experimental 
data are available and for which estimates are made 
in this investigation, is shown in figure 1. Identical 
turret domes are mounted on both upper and lower surfaces 
of a 24-inch-chord section of a symmetrical airfoil 
approximately 19 percent thick. Both turrets are 
rotated -3.19° with respect to the wing chord line and 
the lower turret is located slightly nearer the wing 
leading edge than is the upper turret. These factors 
of turret angle and location necessitate different 
filleting between the dome skirts and the wing surfaces. 
If the flow to which the turrets are subjected is con- 
sidered parallel to the surfaces of the wing, the 
ordinates for the effective turret shapes should con- 
sequently be measured from the wing surface to Include 
fillet thickness and thus effectively change the shapes 
of the turrets by different amounts. 
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Figure 2 shows a comparison of the actual shapes 
of the turrets as installed on the wing, and the effec- 
tive shapes with the wing surface reduced to a straight 
line . 



PROCEDURE 



Pressure distributions for the turrets were esti- 
mated by the methods of reference 1. The boundary layer 
over the streamline shape of the turret was considered 
negligible and the flow was assumed to follow the surface 
without separation. The estimations were limited to the 
three regions considered most important over and near 
the turret in its normal position. In figure 3, these 
regions are identified as: line 1, chordwise over the 
top of and extended in front of the turret central 
meridian; line 2, over the rim or base of the turret; 
line 3, spanwise over the top of and extended from the 
side of the turret at the point of maximum thickness. 

For computational purposes, the effective turret 
shapes (fig. 2(b)) were resolved into simple geometrical 
forms that could be expressed by analytic equations. 
Velocity distributions over the turrets were calculated 
from these equations and interference effects were 
included by adding calculated wing induced velocities to 
the turret velocities; velocities over the wing surface 
near the turret were corrected for interference by 
adding induced velocities due to the turrets. Total 
velocities including interference were then converted 
into pressure coefficients and compressibility adjust- 
ments were applied. 



ESTIMATION OF PRESSURE? OVFR TURRET 



Chordwise central meridian (line 1, fig. 5) .- The 
slope method (reference 1) was used to calculate approxi- 
mate two-dimensional velocity- increment ratios for 
potential flow over simple geometrical shapes best fitting 
the turret profile. By this method,- the velocity dis- 
tribution isa function of the slope with accuracy 
greatest in regions of little slope (near turret maximum 
thickness) and with error most likely in regions of steep 
slope (near turret leading and trailing edges) . 



NACA ACR No. L4G28 



5 



The shapes best representing the turret profile 
consisted of a. circular arc fitting., the forward -portion 
of the profile and a parabolic arc fitting the rearward 
portion with the slopes of the two arcs continuous at 
the point of juncture near the turret maximum thickness. 
(See fig. 2(h).) The origin for equations representing 
these shapes is on the axis of the turret immediately 
below the point of juncture. The velocity-increment 
ratio is expressed as the sum of two integrals corre- 
sponding to the arcs constituting the turret profile in 
the following equation! 



AV 
Vn 



no dy-j/dx -, f T ' E « dy 2 /dx 

I / xT~lE <^ + i dx (1) 

Jl.f.. j o 



where the slopes dy]_/dx and dyg/dx are algebraic 

functions of the abscissa x and AV/V Q is a function 

of the abscissa of any point at which velocities are 
desired x' . In the equation the first integral corre- 
sponds to the circular arc and its limits are from the 
turret leading edge L.E. to the origin or point of 
juncture of the arcs 0. The second integral corre- 
sponds to the parabolic arc and -its limits are from 
the origin 0 to the turret trailing edge T.E. 

Velocity ratios were obtained from the calculated 
velocity- increment ratios by the relation 



X = M + i {2 ) 
v 0 v G + 1 < 2 > 



and are shown as solid lines in figure 4. These calcu- 
lated approximate two-dimensional velocity ratios were 
corrected to the three-dimensional radios ^V/V£j 5. ^ 

f6r the turret by means of' an empirical correction factor 
determined from available data for a similar shape 
(turret B, reference 2) . This factor, expressed as a 
ratio of experimental ( V / V o)3T)' to calculated ( v / v o)g£)j 
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is shown as long-dash lines in figure 4. The magnitude 
of this correction factor was checked for the forward 
portion of the turret by comparing it with ratios of 
calculated velocities for flow over the elliptical 
section of an equivalent oblate spheroid to calculated 
two-dimensional velocities about the corresponding 
elliptic cylinder. The values of ( v / v o).3D obtained 

by multiplying the calculated values of ( v / v o)2D by 

the correction factor at corresponding distances from 
the turret leading edge are shown as short-dash lines 
in figure 4. 

Wing velocity distributions for incompressible 
flow were calculated by the method of reference 3 for 
angles of attack of 0°, 1°, and 4°. Velocity distri- 
butions calculated by this method are usually reasonably 
accurate but experimental values would probably be 
preferable. Figure 5 shows the calculated velocities, 
expressed as interference velocities or velocity-increment 
ratios, for the wing alone at M = 0. 

Total velocities over the turret with wing inter- 
ference were obtained by adding the estimated velocity 
ratios for the turret shape alone (fig. 4) to the calcu- 
lated wing velocity- increment ratios (fig. 5) ; thus 

The resulting total velocities are shown in figure 6. 

Pressure coefficients corresponding to the total 
velocities were obtained from the low-speed relation 

Pi = 1 - (^-) 2 (3) 

and these values were adjusted for Mach number by the 
von Karman relation (reference 4) 
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The estimated peak negative pressure coefficients 
against Mach number are shown in figure 7 with the 
measured values. The intersection • of these curves with 
the curve for critical pressure coefficient P cr (the 

curve of pressure coefficient corresponding to the local 
speed of sound) indicates a critical Mach number of 
about 0.500 for either the upper or lower turret with 
wing interference; this value was taken as the limiting 
Mach number for application of the pressure-distribution 
estimations, since the methods used in the present paper 
do not hold at speeds much higher than the critical. 

Figures 8, 9, and 10 show pressure-distribution 
estimates for the turrets with wing interference at 
angles of attack of 0°, 1°, and 4°, respectively. 
Measured values are shown for comparison. The measured 
values shown in figure 8 were obtained from cross plots 
of measured pressure coefficients against angle of 
attack. 

Rim ( line 2, fig. 3) .- The forward portion of the 
iipper turret can be approximated by an equivalent 
oblate spheroid as shown in figure 11. 

Motion of the oblate spheroid perpendicular to its 
polar axis corresponds to that of the turret. The 
potential is given in terms of the semielliptic coordi- 
nates ti, £, and co in reference 5. The velocity 
relative to the body can be obtained for any point from 




+ 4B- 



1 - u cos to 



(B - B£ cot-H + 0 < 5 > 



t 2 + 1 t 2 + ^ 2 



where 



B = - 



1 
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a is the semimajor axis, and b is the semiminor axis 
for the elliptical section of the oblate spheroid. For 
points on the surface of the spheroid, i = l Q > 

At the surface around the rim (line 2), u. = 0, 
i - t 0 * equation (5) reduces to 

/ v \ = L sln a (7) 
\ °/rim 



where 



and the distance from the leading edge x is given by 

x = a(l - cos oo) 



Over the top (line 1), oo = 0 and 



wh 




x = a(l - \/l - 



The velocity distribution around the rim of the 
forward portion of the turret could be approximated by 
the velocity distribution about the rim of the oblate 
spheroid (equation (7)) as shown by the short-dash lines 
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in figure 12. Experimental data indicate, however, that 
the shape of the velocity distribution around the rim of 
such a regular body bears a close resemblance to that 
over the central meridian in- the direction of motion. 
(See solid line in fig. 12.) At corresponding positions 
along the oblate spheroid and along the turret iength in 
the x-direction, the ratio of velocities at the rim to 
those at the top of the oblate spheroid (see short-dash 
and long-dash lines in fig. 12) were therefore assumed 
to hold for the turret also. Prom this ratio (calculated 
for the oblate spheroid) and the velocities previously 
calculated by the slope method for the chordwise central 
meridian of the turret (see solid line in fig. 12.), the 
turret-rim velocities can then be estimated closely by 
the following relation: 



v 



/_V_\ _ V %blate-spheroid rim /_V_^ 

\° /turret rim fv - ^ \° /turret 

\ %bl 



top 

ate-spheroid top 



Wing interference velocities at a - 1 were added 
to the estimated turret-rim velocities and the resulting 
total velocities were converted to low-speed pressure 
coefficients by equation (3) . Compressibility effects 
were accounted for by the procedure used in the estima- 
tions for the chordwise central meridian. The pressure 
estimates for the rim of the forward portion of the upper 
turret at a = 1° are given in figure 13. Measured rim 
pressures are shown for comparison. 

The rim pressures over the rearward portion of the 
turret could not be calculated. Experimental data 
(reference 6) indicate, however, that for such bodies 
the pressures are nearly constant across the span at 
corresponding distances behind the maximum thickness. 
The pressure coefficients for the rim of the rearward 
portion of the turret were therefore assumed to be the 
same as those along the chordwise central meridian 
(fig. 9(a)) and are shown in figure 13. 

Spanwise at maximum thickness (line 5, fig. 3) . - 
It has already been shown that the forward portion of 
the turret can be approximated by an oblate spheroid. 
(See fig. 11.) 



10 



NACA AGR No. L4G28 



For an oblate spheroid, the velocity over the 
central meridian normal to the flow is equal to the 
constant L of equation (7) . Prom this consideration, 
the velocity over the spanwise section of the turret at 
maximum thickness is assumed to be constant and equal 
to the value previously calculated at maximum thickness 
on line 1. 



ESTIMATION OP PRESSURES OVER WING SURFACE NEAR TURRET 



Velocity decay .- In the determination of loads 
over the wing surface surrounding the turret, the turret 
was considered to be adding interference velocities to 
the velocities for the wing alone. Calculations were 
made along lines extended over the wing surface chordwise 
in front of the turret (line 1 extended) and 3panwise 
from the side of the turret (line 3 extended) as shown 
in figure 3. 

The forward portion of the upper turret was repre- 
sented by the same oblate spheroid used in previous 
calculations (see fig. 11) and velocity ratios were cal- 
culated for this shape by equation (5) with £ = £ 0 « 

Along line 3 extended spanwise, u. = 0 and w = ^; along 

line 1 extended chordwise, u. = co = 0. The calculated 
velocity ratio V/V Q was expressed as the velocity- 
increment ratio AV/V 0 by the relation given in equa- 
tion (2), and the decay ratio AV/(AV) max or 
AV/(AV) min was expressed as the ratio of AV/v Q for 
t, 2 i Q to AV/V 0 for t = S 0 - 

The decay ratio with incompressible flow was calcu- 
lated as a function of t and expressed in terms of 
distance r^ from the turret surface as 

ri = JJEI±.\ (9) 



where a is the length of the semimajor axis. The 
expansion of the field for compressibility was taken 
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into account by plotting AV/(AV) raax against r instead 
of r^ where, from reference 7, : ----- 

. r = ri (10) 
VI - M 2 

Figure 14 shows the velocity decay ratios calculated 
for the oblate spheroid and compared with the values 
derived from measurements on the wing near the surface 
of the turret. The theory of the expansion of the 
field (equation (10)) given in reference 7 is not 
developed for the chordwise direction. The experimental 
values shown in figure 14 were obtained by deducting the 
velocities over the wing alone from the velocities 
obtained at corresponding points on the wing with the 
turret added. 

Spanwise loads (line 5 extended) .- In the deter- 
mination of the estimated spanwise loads, the first 
step was the determination of the velocity ratio for the 
wing alone ^v"/V 0 ) w . Velocity-increment ratios, previ- 
ously calculated for incompressible flow, were obtained 
from figure 5 for the wing alone at angles of attack of 
1° and 4°. These velocity- increment ratios for the 
wing were converted to low-speed pressure coefficients 
by the use of equations (2) and (3) and these values 
were adjusted for compressibility by the von Karman 
relation (equation (4)). The corresponding velocity 
ratios (v/V 0 ^ were obtained from the relation 




(11a) 



The velocity ratio for the turret with wing inter- 
ference (,V/v"o) ma x was determined by use of equa- 
tion (11a) from the estimated pressure coefficients 
over the turret chordwise central meridian (line 1) at 
maximum thickness as shown in figures 9 and 10. 
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The turret interference velocity at the surface of 
the turret was then defined as 




and the turret interference velocity for a point on the 
wing at a distance from the turret was taken as 




where AV/(AV) max is the decay ratio. 

Total wing velocities including turret interference 
velocities were given by the relation 

V _ / V \ , /AV\ 

v o " Ww WA 



Pressure coefficients were obtained from the total 
velocities by use of the formula 




(lib) 



Figure 15 gives estimated pressure coefficients over 
the wing surface spa.nwi.se from the side and over the top 
of the -upper turret for wing angles of attack of 1° and 
4° with .'-tech ir-mbers of 0.275 and 0.500. Experimental 
values are given for comparison. 

Chorlwi se_ loads (l ine 1 extended) .- For spanwise 
loads, ( v / V"o )iy was assumed constant because all calcu- 
lations were made for the same chordwise location. This 
assumption could not be applied to chordwise estimations 
because the wing velocity gradient was significant. It 
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thus became necessary to obtain 



(V/V 0 ) w for each point 



consi'dere'd; otherwise; the procedure was similar to that 
for the spanwise loads. 

The assumption of a stagnation point gave a negative 
turret interference velocity at the surface (turret 
leading edge) 



mated as for the spanwise loads but with the use of the 
chordwise decay ratio given in figure 14(b) and with 
minimum values substituted for maximum. 

The estimated and experimental pressure coefficients 
over the wing surface chordwise in front of the upper 
turret for wing angles of attack of 1° and 4° and at 
Mach numbers of 0.275 and 0.500 are given for comparison 
in figure 16. In the absence of theory concerning the 
effect of compressibility on the extension of the field, 
the decay curve for incompressible flow was used. 



The generally satisfactory agreement shown in fig- 
ures 8 to 10 and 15 between the estimated and the 
experimental pressure distributions may be taken as 
experimental verification of the validity of the method 
of accounting for interference by the addition of 
induced velocities. (This method is of course subject 
to certain limitations that are pointed out in refer- 
ence 1.) In addition, the remarkable agreement between 
the estimated and experimental velocity decay ratios and 
pressures along two lines extending in front of and from 
the side of the turret (figs. 14 to 16) provides an 
excellent check on the potential-flow method of estimating 





DISCUSSION 
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the pressures surrounding a protuberance. As indicated 
by the potential theory, the influence of a protuberance 
on the surrounding pressures decreases very rapidly with 
distance from the surface. 

The effect of compressibility on pressures over the 
turrets is expressed satisfactorily up to the critical 
Mach number by the von Karman relation (equation (4)). 
(See figs. 7 to 10 and 13.) At Mach numbers much higher 
than the critical, however, the von Karman relation fails 
completely to describe the changes that occur. As shown 
in figure' 10, the negative pressure coefficients at 
supercritical speeds may either exceed or fail to reach 
the estimated values. Additional study is needed before 
the pressure loads at supercritical speeds can be pre- 
dicted with confidence. 

The assumption concerning expansion with Mach number 
of the field normal to the turret (equation (10)) is not 
quantitatively checked because the experimental data are 
of insufficient accuracy; from figure 14(a) , however, 
this assumption appears to be at least qualitatively 
correct and the extent of the velocity field outward 
appears to increase at least as much as is given by 
equation (10). It also appears from the experimental 
data (fig. 14(b)) that the assumption of some expansion 
of the field with Mach number chordwise in front of the 
turret would bring the estimation into closer agreement 
with experiment but this condition should not be taken 
as proof that such expansion occurs. 

Although within the range of applicability of the 
method (subcritical speeds) the peak negative pressures 
and consequently the critical speeds have been predicted 
with considerable accuracy as may be seen from figure 7, 
some divergences appear in other parts of the estimations. 
The disparity that exists between estimated and experi- 
mental values in the subcritical range is due largely to 
inaccuracy in calculating the velocities over the turrets 
alone. The divergences for the subcritical Mach number 
range shown in figures 8 to 10, for instance, may be 
explained to a great extent by the fact that the slope 
method used in the estimations is most accurate in 
regions of small slope, as at the maximum thickness of 
the turret, and least accurate where the slope is large, 
as at the front and rear. For the estimations of fig- 
ure 15, constant pressure spanwise across the surface 
of the turret was .assumed as for the case of the oblate 
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spheroid, whereas the actual negative pressure coeffi- 
cients are larger at the center. The agreement between 
estimated and experimental spanwise decay ratios shown 
in figure 14(a) indicates that the departure of the 
shape of the rearward portion of the turret from that 
of the oblate spheroid apparently has very little effect 
on the spanwise decay of velocities. The underestimation 
of pressures over the wing surface spanwise from the side 
of the turret (fig. 15) is accounted for by a slight 



ured pressures chordwise in front of the turret agree 
well although the measured values do not appear altogether 
consistent. 

The general methods employed in the present investi- 
gation may be expected to apply with accuracy sufficient 
for design purposes not only to estimations for similar 
shapes but also to estimations for shapes varying con- 
siderably from the shapes considered herein. Assumptions 
made in connection with the methods, however, should be 
in keeping with the suggestions outlined in reference 1; 
for instance, the assumption of an empirical correction 
factor for reducing calculated two-dimensional velocities 
to three-dimensional velocities may be expected to apply 
to shapes other than those considered herein if the 
ratio of three-dimensional to two-dimensional velocities 
can be determined for similar bodies. A discussion of 
this assumption is given in reference 1. The relation 
assumed between top and rim velocities for the turrets 
in the present report may also be expected to apply with 
reasonable accuracy to other bodies that can be repre- 
sented by simple geometrical shapes, such as oblate or 
prolate spheroids, over which the velocities are calcu- 
lable. The assumption of approximately constant 
velocities spanwise across the rear portion of the 
turrets has been verified experimentally in reference 6 
for other shapes considerably different from those con- 
sidered in the present investigation. 



underestimation of 




CONCLUSIONS 



The applicability of current methods of estimating 
pressures over and near protuberances to turrets located 
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on a wing in regions of high interference velocities 
was investigated. A comparison of estimated and meas- 
ured data resulted in the following conclusions: 

1. Current methods of estimating pressures over 
protuberances were verified. 

2. Methods of predicting the influence of a pro- 
tuberance on pressures over the surrounding surface were 
found to be accurate. 

3. The accuracy of estimations was sufficient for 
design purposes and for predicting the critical speeds 
of the turrets. 

4. The methods used in this paper were not appli- 
cable to estimations at speeds much higher than the 
critical . 
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Figure /.- General details of the turrets mounted on a symmetrica/ wing section. 

All dimensions are given in inches. 
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(a) Actual shapes of the turrets 
as installed on the wing . 



(b) Effective shapes of the turrets 
with wing surface reduced to 
a straight line . 
figure 2.- Actual and effective turret- profile shapes. 
All dimensions are given in inches. 
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Fig. 4a, b 
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Figure 4-.- £3 f /mated vs/oc/t/es far the chord wise 
central meridian of the turrets without wing 
interference. /W- o . 
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figure S.- Interference ve/ocities 
calculated -for the winq a/one . 
A/7 — O . 
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Fig. 6a, b 
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Figure S. - Estimated velocities -for* -the 
chordw/se centra/ meridian of the turrets 
with wing interference. M = O . 
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Fig. 7a 
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Figure 7. - Variation or" peak negative pressure 
coefficients with Mach number for the 
chord wise centra i meridian of the turrets 
with wing interference . 
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(a) Upper turret. 
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Figure S. - Rressurs distributions for the 
chorduvise centra/ meridian of the turrets 
with wing interference. 4C - o° . 
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Fig. 9a, b 
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. I L I i J I i 1 ! 1 1 

O I Z 3 4 S 6 

Distance from turret ieqding edge, in- 



Figure S.~ Pressure c/i'sfributions for the 
chordwise centra/ meridian of f/iB turrets 
with wing interference. <X = l " • 
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Fig. 10a, b 
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Figure iO.~ Pressure distributions for the 
chord** is e centra/ meridian of the 
turrets Hsi+h wing interference . c(.= 4- a . 
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Figure II." The forward portion of the upper turret represented by an oblate spheroid. 

All dimensions are given in inches . 
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figure 12. . - Estimated velocities around the rim 
or the upper turret without wing interference. 
Avf - O . 
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Fig. 13 
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Figure /3.- Pressure distributions for the 
rim of the upper turret with wing 
interference . ec = / • . 
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Fig. 14a, b 
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(a) Spanwise from the 
■\ , , Side of the turret. 
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Figure 14-.- Velocity decau with distance from 
the turret- surface . 



NACA ACR No. L4G28 



Fig. 15a, b 
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figure /5.- Pressure distributions Spanwise over 
the upper turret at maximum thickness and over 
the wing Surface near the turret. 
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Figure 16.- Pressure distributions chordwise over the wing Surface 

in front of the upper turret . 
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